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Abstract 

We study the dynamics of a single hole in an otherwise half-filled two- 
dimensional Hubbard model by introducing a nonlocal Bogolyubov transfor- 
mation in the antiferromagnetic state. This allows us to rewrite the Hamil- 
tonian in a form that makes a separation between high-energy processes 
(involving double-occupancy) and low-energy physics possible. A diagram- 
matic scheme is developped that allows for a systematic study of the different 
processes delocalizing a carrier in the antiferromagnetic state. In particular, 
the so-called Trugman process, important if transverse spin fluctuations are 
neglected, is studied and is shown to be dominated by the leading vertex cor- 
rections. We analyze the dynamics of a single hole both in the Ising limit and 
with spin fluctuations. The results are compared with previous theories as 
well as with recent exact small-cluster calculations, and we fing good agree- 
ment. The formalism establishes a link between weak and strong coupling 

methodologies. 
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I. INTRODUCTION 



Strongly correlated fermion systems have been the focus of interest in the last few years 
following the discovery of high temperature superconductors. At the heart of the problem 
is the interplay between the itinerant character of the charge carriers and the existence of 
localized spins (on some time scale) imposed by the single-occupancy constraint. Consid- 
erable effort have been invested in the understanding of the Hubbard model and its strong 
coupling cousin, the so-called t-J model, in particular in two dimensions. At half filling, 
the strongly correlated Hubbard model reduces to the Heisenberg model, which is now be- 
lieved to exhibit antiferromagnetic order on bipartite lattices in two dimensionsB. Further, 
for arbitrary correlation strength U, the Hartree-Fock solution leads to a commensurate 
spin-density wave (SDW) state which approaches the Neel state as U increased. In addi- 
tion, a straightforward RPA calculation for quantum spin fluctuations has produced a spin 
wave mode which in the large U limit coincides with what is expected from the Heisenberg 
modeil. 

Away from half filling, no consensus exists yet except the Nagaoka theorem. In the 
weak-correlation limit it is known that doping tends to induce an incommensurate SDwi. 
In particular, the antiferromagnetic state starts to be discommensurated at small doping by 
deforming into a domain wall structureB^i. Increasing the strength of correlations suppresses 
the amplitude fluctuations, enhancing instead the instability in the transverse channefl In 
particular, a mean field analysis for the t-J model has predicted a spiral phase0. However, 
inhomogenous phases have also been argued to be stable at least for small doping even in 
the large-U limiti'i. So far few substantial studies have been devoted to exploit the effects 
of quantum fluctuations or to investigate in detail the crossover between weak and strong 
correlation in the doped case. 
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Given the different methodologies used in discussing the weakly and strongly correlated 
regimes of the Hubbard model, it is not at all obvious however why the smooth interpolation 
should occur. In particular, since weak coupling methods are based on the assumption that 
the coupling is much less than the bandwidth 8t, it is unclear how the local constraint of no 
double occupancy is automatically respected in extrapolating the weak coupling methods 
to the strong coupling regime. The purpose of the present work is partly to fill that gap. To 
moderate our ambition we shall limit ourselves only to the case of a single carrier. This will 
allow us to assume the two-sublattice structure in our analysis, since the antiferromagnetic 
order is then believed to persist. Moreover, currently available data on finite clusters make 
reliable quantitative comparisons possible. 

The calculations are based on the following idea. By introducing a nonlocal Bogolyubov 
rotation on the Hubbard Hamiltonian we set up a new basis in which we can readily see the 
differences and connections between weak and strong coupling approaches. In particular, 
the low energy scattering processes can be diagrammatically isolated from those at high 
energy without ambiguity. This in turn makes it possible to show how the latter processes 
are scaled away as U/t becomes sufficiently large, thus virtually excluding the possibility 
of double occupancy in the strong coupling limit. 

The paper is organized as follows: in section |l| we introduce the Bogolyubov rotation 
and derive the effective Hamiltonian within the new basis. In section |ITT| we take as example 
the half filled case to examine quantum spin fluctuations. The analysis will be performed 
entirely within the rotated basis in which (1) the Hartree-Fock ground state has exactly 
the same spin configuration as the unperturbed "noninteracting" state, namely the SDW 
state; (2) fluctuations above this SDW state can be classified and diagrammtically distigu- 
ished according whether double occupancies are involved or not. Such a distinction from 
the conventional RPA analysilEI proves to be an important advantage. Based on the 
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observations made in section |T|, the Feynman diagrammatic scheme is then apphed to the 
problem of a single hole in the strongly correlated regime in section |V|. We shall discuss in 
detail the effects of fluctuations on the dynamics of a single hole and compare with results 
derived using other available techniques. Comments and a discussion are given in the final 
section. 



II. THE BOGOLYUBOV TRANSFORMED HAMILTONIAN 

The standard Hamiltonian of the Hubbard model is 

H = -tJ2 + ^-c-) + ^E^^'T^^i ' (2.1) 

where the Ci^^ are annihilation operators for a fermion at site i with spin a, Ui^a is the 
corresponding number operator, and j) indicates summation over the nearest neightbor 
bonds of a square lattice, each bond being counted once. The ratio between the onsite 
repulsion parameter V and the hopping energy t determines, together with the bandfilling, 
the physics of this model. In this section we introduce the Bogolyubov rotation of the 
fermion operators and then derive an effective Hamiltonian for the strongly correlated case 
{U ^ t) near half filling. 

We begin by going to momentum space and rewriting the Hamiltonian within the re- 
duced Brillouin zone (RBZ) defined by + |^y| < tt. This can be done by introducing the 
spinor notion "^^ik) = {Cl^, C^^q^), where Q = (±7r, ivr). We then obtain (see Appendix 
A) 

ka 

kkiq 

+ [^{k + q)a^^>^{km\{k' -q)a'm^{k')]} , (2.2) 
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where the a* {i = 1,2,3) are the standard Pauh matrices, and = — 2t(cos fc^,' + cos ky). 
One should notice that while the summation over momenta runs only over the RBZ, k±q 
can freely go beyond the boundaries of the RBZ. Each state is then easily seen to be counted 
once and only once. 

Next we set up a new basis by introducing the following unitary transformation of the 
fermion operators: 



( 



Uk Vk 

and for Uk, Vk we take the large-f/ expansion 
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(2.3) 



This of course is the Bogolyubov transformation used in the weak coupling approach in 
diagonalizing the Hartree-Fock Hamiltoniani. However, we emphasize that the operation 
here will be applied to the full Hamiltonian instead of the mean field one, and therefore 
is exact irrespective of the form of the coefficients Uk-, Vk- As it will soon be clear this 
procedure will permit us to examine to what extent the weak coupling approach overlaps 
with the strong coupling methods. 

The transformed Hamiltonian is rather involved when written in momentum space and 
includes both interband and intraband scatterings as well as conditional single particle 
interband hoppings. To render things more transparant we turn instead to real space. 
Assuming U/t^l and keeping only terms up to 0{t'^/U^), we can easily rewrite eq.( p.3| ) 
in real space and obtain a nonlocal transformation: for the site i on the even sublattice and 
cr =1 or site i on the odd sublattice and cr =| one has 

^io- = ''^i ~ Tj'^ ~ ^ '^i+S-5' 5 (2-5) 

5 S,5' 



whereas in the opposite case 

-'^+rj^ - ^ E '^+s-S' ■ (2.6) 

s S,S' 

and J = At'^/U as usual. The main difference between the new basis and the original one 
is thus very clear: the electrons are now classified in terms of the category (c or v) they 
belong to, instead of spin: a v electron has spin J, if it sits on the even sublattice and spin 
t if it sits on the odd sublattice, and vice versa for the c-electrons. Thus, a fully occupied 
v-band is essentially a Neel ordering in the original basis (namely, we observe a downspin 
on the even sublattice and an upspin on the odd sublattice). We note that the particles 
created by the r-operators are not strictly localized on their respective sites, but are rather 
in the Wannier states corresponding to the mean-field Neel state and thus are to a certain 
degree delocalized on nearest (amplitude oc t/U) and next-nearest (amplitude oc {t/UY) 
neighbors. 

The effective Hamiltonian is finally obtained by substituting the above tranformation 
into the Hubbard Hamiltonian. We thus have, retaining only terms up to 0{t^ /U): 

H^H, + H2 + Hd , (2.7) 
+ 2J5]« - <) - jJ^iS^S] - \n,nj) + ^ J^i^jS- + S'S}) , (2.8) 

- I E i^f^'^+srt'rl + h.c.) , (2.9) 

i,S^S' 
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i i,55' 

- I + h.c.) + t J2Krf^1+sri^+s + h.c.) . (2.10) 

i,5S' iS 

The spin index in the original Hamiltonian now is replaced by the v and c band indices, as 
discussed following eqs. ( |2.5| ) and (|2.tjD. The spin operators in Hi are defined by Sj = rfr^ 
on the even sublattice and S~ = rYrJ on the odd sublattice and represent local spinfiips 
in the new basis. Further, fiifij/A — S^Sj in Hi is a rewriting of (n^n^ + n^nJ)/2 which 
represents nearest neighbor intraband scattering. We have expressed the Hamiltonian in 
three terms, namely Hi, H2, H^-, to emphasize the fact that all interaction processes in 
Hd involve double occupancy, and in H2 the electrons will always go beyond the nearest 
neighbors. We also notice that in Hi and Hd there are terms involving six fermionic 
operators. These terms in fact cancel (not surprisingly, given that the original Hamiltonian 
only involves terms containing up to four fermion operators and that the transformation to 
the new operators is linear), however, we find the forms ( p.8| ), ( p.lO| ) rather useful because 



they clearly separate terms involving double occupancy from those which don't. 

It is clear that at half-filling (one electron per site) the "non-interacting" (or Hartree- 
Fock) ground state of the tranformed Hamiltonian (|2.7] ) is not the paramagnetic Fermi sea 
as in the original basis, but rather a state with all the w-states filled and all the c-states 
are empty. That is to say that in the new basis electrons will first fill the f-band in which 
the spin of each electron is seen to point down on the even sublattice and up on the odd 
sublattice. The Hartree-Fock energy of this state is Eq = —NJ. 

The Hartree-Fock quasiparticle energies are obtained by calculating the average value 
of in a state with either a c-particle added or a f-particle taken out. One then finds 

Hhf = 5^(efcr^V^ + arl^rl) - NJ , (2.11) 
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where = — ^ + f^l/f^ and Ck = J — ^l/f^, and the fc-space operators are defined by 

r^''' = ^yr,^'V*=-^ , (2.12) 

with the fc-summation covering the full Brillouin zone. We note that the operators thus 
defined are not eigenfunctions of cr^, because the spin projection of e.g. a f-particle is 
different on the even and odd sublattice. Spin eigenstates can however be straightforwardly 
constructed noting that C.'j^^ = ^fc+g- The operators with well-defined spin then are r^'^ = 
(r^'^' ± r'^^o) / , where A;-space is now limited to the reduced Brillouin zone. In terms of 
these operators the Hartree-Fock Hamiltonian takes the form 

Hhf = + arllrL) - NJ . (2.13) 

ka 

The c-band is now seen to be separated from the fully occupied f-band by a gap of order 
U. In addition, the Hartree-Fock ground state at half-filling has a fully filled f-band and 
an empty c-band corresponding to Neel spin ordering. 

Going back to the full Hamiltonian ( p^.T] ) it is clear that for electron density n < 1 the 
c-band is drawn in only through interband interactions. To get an idea how the c-band 
becomes involved, we notice that interband interactions can actually be divided into two 
classes. Interband processes described in can bring in the c-band but always provoke 
double occupancy simultaneously. At the Hartree-Fock level this of course pushes the c- 
band by U above the f-band, and it is clear that beyond Hartree-Fock, there will always 
an intermediate state of energy ^ U involved. On the other hand, interband processes 
contained in Hi + H2 can actually transfer the f-band electrons to the c-band or vice 
versa at much less energy expense. One such example is the interchange of two nearest 
neighboring f-band particles. This amounts to having a spinfiip at each site since these 
two particles have opposite spins, or equivalently the pair of f-band particles hops to the 
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c-band in the end. Such a process is reahzed via the interband pair hoping described by 
the last term in Hi and costs only an energy of order 0{J). 

It is physically clear that scatterings involving double occupancy will actually be cir- 
cumvented in the strong coupling limit and therefore render no substantial contribution to 
the low energy physics. Nevertheless, this is not so obvious if we deal directly with the 
Hubbard model using Feynman's diagrammatic technique. It turns out and will be shown 
below that the advantage of the Bogolyubov rotation proposed here is that it provides a 
basis in which we are able to clearly demonstrate with Feynman diagrams the irrelevence 
of Hd in the effective low energy properties. 

Once established that will eventually be scaled away, we are left with Hi + H2 to 
describe low energy physics of the large-C/ Hubbard model. In particular. Hi is scarcely 
different from the familiar t — J model, except for the interactions in the (charge and 
spin) density channel. Here we note that the "noninteracting" state we start with in the 
new basis exhibits Neel ordering. The spin fluctuations above the Neel state are then of 
ferromagnetic nature in the longitudinal spin channel. On the other hand, since the spin 
rotational symmetry has been retained in the transformation fluctuations in the transverse 
spin channel are described in exactly the same manner as in the t-J model. 

III. QUANTUM FLUCTUATIONS AT HALF FILLING 

We devote this section to the discussion of quantum spin fluctuations in the half-flUed 
Hubbard model. The flnal result is actually known either from spin wave calculations or 
from direct extrapolation of weak-coupling results. The convergence of these apparently 
unrelated approaches raises the question of a possible smooth crossover between the two 
extreme coupling regimes. The main point which remains unclear is why and how the 
local constraint of no double occupancy becomes automatically respected when one goes to 
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large U within the mean-field formalism which is frequently believed to be limited to weak 
coupling. The purpose of our presentation is to demonstrate that within the new basis 
we are able to show explicitly and in Feynman diagrams how the higher energy processes 
involving double occupancies {i.e. Hd) are in the end scaled away and become irrelevant to 
the spin dynamics. 

The antiferromagnetic ordering implies the the existence of gapless spin wave excitations 
in the transverse spin channel. This can be analyzed by studying the dynamical spin 
susceptibility in the presence of Neel ordering which is now defined as a 2 x 2 matrix: 



X"+(x,x',t) 



^ {T (5- (t)^+(0))) (T (^- (t)^+ (0))) ^ 



(3.1) 



As convention, the position coordinate x is denoted by Xe if the spin is on the even sublattice 
and by Xq if on the odd sublattice. We have added a bar over the spin operators in order to 
emphasize that we are probing correlations between spins defined in terms of the original 
C-operators (as opposed to the ^-operators occurring in eq.(p.8|)). 

The calculation for the above matrix at the Hartree-Fock level is straightforward. We 
first transform the susceptibility matrix such that S^, are written in terms of r-operators 
as presented in Appendix B. The leading contribution is then readily obtained. Up to 
0{J/U) only the diagonal terms in eq. (|3.1j ) are nonzero and are given by 

XHFuiQ, ^) = XHF22 = ]^ I] ^ _ ^ ^ " (3-2) 

With little algebra0 one can verify that this agrees with the result^ of SWZ up to 0{J/U). 

To calculate the dynamical susceptibility in the presence of interactions we need to 
take some care. In all previous studieJiE, the RPA scheme was adopted within the 
paramagnetic C-operator basis in which the interaction in the transverse spin channel is 
exactly U. Here we rather write the original spin operators (S) in terms of the r-operators. 
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Then the simplest Hartree-Fock ground state aheady has the correct antiferromagnetic 
structure, and the interaction between the Hartree-Fock particles is directly obtained from 
eqs.( |2.7| )-( p.lOD . To facilitate the following discussion we write the interaction explicitely 
in momentum space. The terms relevant for the transverse spin fluctuations are 



kklq 

+V i — > c) , 



(3.3) 



with 



gi{k,k') = U 



ei + e 



k+q 



+ el + el_ 



2f/2 



u 



(3.4) 



(3.5) 



We note that no 0{t) terms enter as can easily be understood for half filling. Neither 
do the intraband interactions enter since they contribute only to longitudinal fluctuations. 

is now composed of two parts. The gi part is equivalent to the first two terms in Hd 
and describes interband scatterings involving double occupancy. On the other hand, the g2 
part describes the transfer of a pair of particles between two bands and corrsponds to the 
last terms in Hi and H2. This part will be shown to be the one responsible for the spin 
dynamics. Finally it will be helpful to note that up to 0{J/U), gi and g2 can actually be 
factorized by gi = Uh{k)h{k')lA and g2 = U{h{k)f2{k') + /i(fc')/2(fc))/4, with h and 
being defined in appendix B. 

We are now in the position to study the transverse spin fluctuations in the SDW state. 
Following appendix B we shall calculate 



2N 



Y.fi{k)h{k') 



kk' 



( \ 

ail cti2 
0.21 ^22 



(3.6) 
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to get the dominant contribution. The definition for the matrix elements is given in ap- 
pendix B. A shght comphcation arises in the calculation due to the fact that in the inter- 
action Hamiltonian two interaction vertices gi and g2 are involved, and both of them 
are momentum-dependent. However, we should notice that the gi interaction conserves 
the number of particles in each band. In other words, for an excited particle-hole pair 
of opposite spin the ^fi-term contributes to multiple scatterings between the particle and 
hole. On the contrary (72^interactions give rise to destruction and formation of p-h pairs of 
different spins. Within the RPA scheme we then find that gi and (72 will operate distinctly 
and the two interactions do not mix up in Feynman diagrams. A typical diagram for 
incorporating both interactions is shown in fig.|l|. 

Having learnt that the gi- and (72^type interactions do not mix up at RPA level, we can 
deal with them separately. Let us look at the (yfi-interaction first. As is evident, the off- 
diagonal terms in eq. ( |3.6| ) remains zero with only gi interactions. To get the diagonal terms 
of eq. (|3l6|) , we need to evaluate the diagram illustrated in fig.|l|b. After a straightforward 
calculation we obtain 



X-,,\,M.uj) = . , (3.7) 



where 1/A{k, q) = ±uj — C.'f.^g + 'Cfc + '^0^- particular, taking the limit of U ^ t gives the 
(yfi-renormalized susceptibility as 



y x'iq, J 



(3.1 



We notice that going from eq. (|3.2|) to eqs. (|3.8|) , ( p.9|) the high energy scale of 0{U) has 
disappeared due to the renormalizaion by multiple gi scatterings and has been replaced by 
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2J which is the energy one needs to create a local spin flip. This is natural because after 
all a local spin flip can occur without involving double occupancies. 

Finally we switch on the g2 interactions. As explained before and shown in flg.|l| within 
RPA they have no effect on the structure within the bubble but play their role by forming 
junctions between (yfi-renormalized bubbles. Besides, g2 interations also render a nonzero 
contribution to off-diagonal terms and 021 in eq.( p.6|) , as is shown diagrammatically 
in flg.|^. Again the work involved in the calculation is routine since g2 also factorizes as 
(72 (A;, k') ~ (/i(fc)/2(fc') + /i(fc')/2(fc)). We only mention that because of the way in which 
g2 is factorized, we shall meet in the intermediate states bubble diagrams like flg.^D but 
with vertex parts /i or /2 attached to one of its two legs. For example, when fi is attached 
to one leg and f2 to the other, the diagram gives ^J'^q/ {U{±uj — 2J)), whereas /2 is attached 
to both legs leads to 4J^7q/(t/^(±c<j — 2J)). The flnal result for the susceptibility then is 
obtained as 



( . . . \ 



1 + — 

^ 2J 



\ ^ 2.J ) 



(3.10) 



where ujq = 2Jy/l — is the standard spin wave spectrum. The same expressions can 
be found from Singh and Tesanovic's resultsUll by directly extrapolating the RPA results. 
However, with the help of Bogolyubov rotation and by working within the rotated basis we 
have explicitly shown here how the higher energy scale gets irrelevant in the spin dynamics 
within the framework of a diagrammatic approach. 

At last, we point out that the above results can equally be derived by considering only 
Hi and neglecting H2 and entirely, using the same Feynman diagrammatical scheme. 
From the Hartree-Fock solution of Hi we immediately get eq.( ^.8D , and the flnal result 
of eq.( 3.1(]|) is obtained following RPA for the Heisenberg exchange interactions in the 



transverse spin channel. The coincidence teaches us two things. First, it shows that the 
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t — J model represents the Hubbard model more faithfully in the strongly correlated regime 
than the t — J — J' or t — t' — J models. Secondly, one sees that it is possible to extend the 
Feynman diagrammatic technique to the strong coupling regime by working directly within 
the sub-Hilbert space of no double occupancy. In the next two sections we shall use this 
fact to discuss the dynamics of a single hole in the otherwise quantum antiferromagnetic 
background. 

IV. THE DYNAMICS OF A SINGLE HOLE 

We now begin to examine the motion of a hole in the antiferromagnetic environment. 
For this purpose we take away one downspin electron from the v-band and study the 
Green's function for the t'-band hole 

g:{k,u;) = (Okfel ^_ff\^o+ ^L|0) (4.1) 

within the sub-Hilbert space of no double occupancy, as defined in the last section. Here |0) 
represents the full v-hand, i.e. the Neel state. For the reason given at the end of the last 
section, and also to render possible a direct comparison with available numerical results 
obtained notably for the t — J{Jz) model we shall in the following (and unless specified 
otherwise) neglect H2 which always leads to jumps beyond the nearest neighbors. The 
principal subject we shall deal with is therefore Hi. 

To begin we shall also neglect spin fluctuations altogether {t — Jz model) . The Hartree- 
Fock state of Hi is defined by two degenerate energy levels {u^ — and Uc — '2,J respec- 
tively). By ignoring spin fluctuations the only remaining interaction in Hi is the 0{t) 
hopping term which takes effect when a hole is introduced and is allowed to move: 

Hi = - <)^'V+5(i - <+5) + (1 - <+.)rrVr(i - <)} . (4.2) 
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Obviously, this term correlates the hopping motion of a particle (or hole) with the occupa- 
tion of its surrounding sites and thus is responsible for the coupling between the moving 
hole and its Neel environment. We now notice that Hj only destroys particles on singly 
occupied sites (double occupancy of a site means that there are both a c- and a f-particle) 
and only creates particles on empty sites. Therefore acting on a state containig only empty 
and singly occupied sites as in ( [4.1|) one can rewrite Hj in the simpler form 



Hi = - O^f^-^s + <V?(1 - <)} ' (4-3) 

which now contains four-fermion terms only. For later convenience we write down its 
Fourier transform in the RBZ explicitely: 

= -| Y.^^^<^l-.i^t'-qfl', + ^^-iVl'+.T-ll-l-.i + h.c.} , (4.4) 

kk'q 

where = —zt'jk = — 2t(cos(fca;) + cos(fcy)). We note in particular that the momentum in 
Ek is always that of the the f-band hole which has the same spin as the c-particle. This 
fact will be used on numerous occasions when we need to determine vertex functions. In 
addition, the hole notation f = r"!" has been introduced in order to facilitate the discussion. 



A. The incoherent spectrum 

It is evident that the first nonzero contribution arises only in the second order correction 
to the selfenergy. Naively one would have expected, as is readily derived from the interaction 
Hamiltonian, contributions from the two diagrams shown in fig.|^. Nevertheless, in the 
absence of spin fluctuations the contribution from fig.^D vanishes since a the lower {v, t) line 
represents hole propagation, which implies that the (c, |) line should represent creation of a 
hole in the c band. This of course is impossible, the c band being initially empty. In other 
words, a pair of a c-particle and a v-hole of opposite spins can not appear before the hole 
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starts moving. Consequently for an initial downspin hole we only see in the intermediate 
state an upspin hole, accompanied by a p-h pair formed by c-band upspin particle and 
f-band downspin hole. The diagram is then evaluated according to the standard Feynman 
rules and is given by: 

^ik,co + ^0^)=^^J2^l^ ^g;^'\k,,u;, + ^0^)x^k-k^,co-co,) . (4.5) 

fcl 

The hole notation has been emphasized in the expression hj u + iO~^ such that the Green 
function of a hole below the Fermi sea is seen as a retarded function of time. The noninter- 
acting Green function is therefore Q""^^^ = l/(uj + iO^), the p-h bubble is easily determined 
and is given by eq. (|3.9|) which is also a retarded function. Moreover, the vertex functions 
at two intertacting lines are all determined by the momentum ki of the scattered hole in 
the intermediate state. We thus obtain as final result: 

The factor z reflects that the hole moves to its z nearest neigbors. The absence of momentum 
dependence corresponds to the fact that the approximation actually binds holes to its origin 
by allowing only back and forth motions towards its nearest neighbors (as implied by the 
momentum summation over 7^). Along this line, the p-h pair in fig.^ is easily seen to 
represent the local spinfiip the hole creates at its initial position once it hops. 

A more complete picture arrises when we introduce selfconsistency in fig.^. The proce- 
dure is immediately seen to permit the intermediate upspin hole to move further forward. 
The new Feynman diagram (fig.^ requires us to solve the following selfconsistency equation: 

fcl ^ 

After substituting the bare p-h bubble, eq.( p.9|) , the coupled equations for S and Q are 
further reduced to a simple iterative equation for the spin and momentum independent 
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Green function: 



= ^ • ^4-8) 



This agrees with the resuh of Shraiman and Siggia for the incoherent hole spectrunJi^ii^ 

The physics is rather clearly reflected in the diagrammatic picture. Because of the 
selfconsistency, the selfenergy is now described by a series of Feynman diagrams of order 
2,4,6,.... At a typical order 2n (fig.^), a sequence of n bubbles is first formed and the 
bubbles are then eliminated one after another in reversed order. Since a momentum- 
independent bubble of opposite spins in the Feynman diagrams means a local spinflip 
in real space, the above sequence describes actually a string of overturned spins in the 
intermediate state. In addition, the number of bubbles n is the length of the string for that 
particular order. In conclusion the selfconsistency has, on the one hand, made possible 
to allow for infinitely many particle-hole excitations in the intermediate states along the 
way the hole hops forward; on the other hand it guarantees that the hole can only follow 
the self-retracable paths defined by Brinkman and Ricelii and remains bound to its original 
position. We have thus properly incorporated the contributions from incoherent excitations 
within the diagrammatic scheme. 

The numerical solution of eq ([4.8|) defines a ladder spectrum or a sequence of bound 
states, with the lowest state being order of t{J /t)"^/^ above the lower edge of the Brinkman- 



Rice band cuq = —2^t. The discrepancy with the original — It of Brinkman and 

Ricelii is understood according to Kane et alEi: in summing 7I in eq. (|4.71 ) contributions 
from all the nearest neighbors are assumed at every stage of the intermediate states. This 
results in overcounting from the second step on, since in fact only forward going steps should 
be counted0. Finally we mention that inclusion of H2 in the above calculation would have 
destroyed the ladder structure because then the hole is already mobile at the Hartree-Fock 
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level. Strictly speaking, this term is of course present (and of order t'^/U) in the Hubbard 
model (but not in the t — J model). 



B. The coherent spectrum: Trugman's motion 

So far the hole is completely localized. In the absence of spin fluctuations it becomes 
mobile only when it manages to move without disturbing the antiferromagnetic background. 
As Trugman has observed^, the simplest path to permit this is one where the hole hops 
around a plaquette one and a half times. The hole then moves to its nearest neighbor 
on the same sublattice while cleaning up the string of reversed spins in the end. In the 
diagrammatic scheme, we need to resort to vertex corrections to account for such paths. In 
fact, earlier studies have realized that Trugman's path corresponds to the leading nonzero 
vertex correction involving two loopsBEl. Within the present framework, the same process 
is described by a unique diagram as shown in fig.^. The calculation for this bare 6th 
order vertex correction is straightforward noting that both the hole lines and the particle- 
hole bubbles in the diagram are momentum independent and retarded in time. A slight 
complication occurs only due to the momentum dependence in the interaction vertex. For 
our purposes we only need to keep in mind that the vertex function at each interaction line 
in the diagram is determined by the momentum of the scattered hole in the intermediate 
state (the bold line in fig.§, as already encountered in obtaining eq.( [4.5|) . The rules for 
calculating the diagrams are again the standard Feynman rules. The result for the bare 
6th order correction is then found to be 



q 



The denominator in the above expression can easily be explained following fig.p3. The 
hole loses an energy of 2 J each time a p-h pair or a spinfiip is excited, its energy reaches 
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— 6 J after three emissions and is only regained via three more steps in which it absorbs 
the emitted pairs (the spinflips are immobile once created). The maximal crossing of loops 
explains the order the hole must follow in the last three steps, the only possible path at 
6th order is then the well-known Trugman loop (fig.^). The delocalization of the hole is 
shown by the momentum dependence of the renormalized selfenergy ( ^.91 ). 

Obviously with closed paths of longer lengths the hole is able to move further away from 
its initial position. Evidently for any translation not affecting the spin structure it can only 
end up on the same sublattice, which in turn indicates that only vertex corrections of order 
4n + 2 will contribute. The next nonzero contribution then appears naturally at 10th order. 
Not supprisingly, the number of possible paths grows considerably as one goes to higher 
orders, fig.0 describes the three processes responsible for the hole delocalization at 10th 
order. Apparently, the first two processes share a common feature in that in both cases 
the hole attempts to create a maximum number (here 5) of spinflips before it eliminates 
them again. They differ from one another only by the fact that they follow distinct paths 
in eliminating the spinflips. fig.|^ has maximal crossing of loops and is obviously the direct 
extension of the 6th order correction. Its real space configuration is therefore described 
by fig.^. The hole follows a fixed path {i.e., 1-2-3-4-5-6) and hops two loops minus two 
steps to its diagonal neighbor (i.e. its first neighbor on the same sublattice) while erasing 
the trace left behind. On the other hand, in fig.0b the hole changes its path starting from 
the 8th step and ends up at one of its second nearest neighbors on the same sublattice (see 
fig.^D). Finally in the third process (fig]^c) we see at most three spinflips in the intermediate 
states and the hole runs along the path illustrated in fig.|]c and also moves to its second 
nearest neighbor sites. 

To get an idea of how those higher order corrections affect the hole mobility, at present 
we will limit ourselves to the first category of closed paths. In other words, we shall only sum 
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the series of maximally-crossed diagrams of 4n + 2 orders (n = 1,2, ...) as demonstrated 
in fig.^ to get a picture of the mobile hole under Trugman's mechanismEl. In a short 
while we shall actually show that as far as low energy physics is concerned, the dominant 
contribution comes only from the leading (6th) order vertex correction. Fig.^ is evaluated 
along the same lines that lead to eg. (^.91) and we get: 

E^k, u + zO+) = ^Y1 ^(9' ^ - 2^)^9K + ^(^' 9)} • (4-10) 

We have introduced selfconsistency in order to simultaneously account for the self-retracing 
paths discussed in the last section. Thus, in the absence of V{k, q) this selfenergy just re- 
produces the incoherent result, eq. ([4.8|) . V{k, q) sums the vertex corrections in fig.^. Using 
the bare Green functions which are momentum-independent the momentum dependence of 
V is entirely due to the momentum dependence of the factors Sk in eq.( [4.4|) . We then have 

V{t^.q) - 2^ 2J){uj - 4J)2(cu - 6J) 

92 

N ^ {Lj-2j){u-4jy{Lj-6jy{Lj-8jy{u-ioj) ^ ■ ^ ■ ) 

92 

To help to understand this equation we note that ^2 the above expression is the hole mo- 
mentum just one step before a complete closed path (see the diagonal hole line in fig.^. In 
addition, the only momentum dependence comes from the interaction vertex, and summa- 
tions over all but one internal momenta have been performed noting that —tSk is obtained 
by summing q in SqEk-q- The omitted higher order terms in ( [4.11|) can be obtained straight- 



forwardly remembering that hare hole lines and bubbles are used inside the vertex part. The 
final result is determined by solving the following equation: 

V{k^ = E 2i^^F(fc', , (4.12) 

fc' 
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Without going into the details of a numerical solution, we can immediately learn two 
points by examining the above equation. First, as long as J/t ^ 1 close to the bottom 
of the Brinkman-Rice spectrum (tu ~ cjq = — we have F{k,uj) ^ {zt'^'-fkY /uj'^. In 
other words, the higher order vertex corrections are quantitatively insignificant around the 
minimum of the spectrum since every additional two loops add a factor of roughly 1/16. 
Therefore, the mobility of the hole is dominantly determined by the leading vertex correc- 
tion. The second point needs some elaboration. The statement is that in any event, the 
vertex correction above alone is not expected to give quantitative satisfaction. Since the 
point has in fact already shown up in eq. ([4.9|) , we take the case of the leading vertex cor- 
rection to demonstrate the essential fact. Recall that in obtaining the preceeding equations 
it is always implicitly assumed (by using eq.(|3l9|)) that the energy cost to create a spinflip 
is 2 J. This is not correct in the presence of a hole. From fig.|^c we actually see that in fact 
it costs the hole only 3 J/2 to perform the first step, J for the next, and J/2 for the third 
step, and vice versa for the energy-recovering steps. The evolution of the energy cost for a 
hole moving around the plaquette is shown in fig.p!0|a. In the same manner, the J factors 
in both eq. ( |4.11| ) and eq.( [4.13| ) are seen not to give a correct description of the energy cost 
in the intermediate states. As can readily be understood this failure is directly related to 
the mean field approximation we have applied in the longitudinal channel (the third term 
in Hi): the intraband correlations between holes at different intermediate states have so 
far been neglected. The renormalizations due to the correlation effects in this channel can 
be easily taken into account at fixed order in the vertex corrections, as done in eq.( [4.14|) 
below. However, a systematic procedure to all orders seems to be difficult to find because 
the energy expense in the intermediate states of a closed path is configuration-dependent. 
Take the 10th order maximally crossed diagram as example: we see the energy evolution 
shown in fig.p!0|b which is apparantly unrelated to what we have observed in fig.p!0|a. This 
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situation is in sharp contrast with the Brinkman-Rice case where we know that except for 
the first step the energy cost to extend the path by one step is always J. 

After these remarks we now concentrate on the leading order vertex correction. This 
can easily be done by truncating the iteration eq.( |4.13| ) by its first order approximation. 
In order to take into account the correlations in the longitudinal channel we also include 
renormaliztions so that the J factors agree with the analysis of the last paragraph (fig.|lO|a). 
Formally this can be achieved (at fixed order in Hj) by using the second and the third term 
in Hi as the zeroth order Hamiltonian. The final equation we then have is 

S^(fc,^+.0+) = ^E^(^'^-3^/2){7?+]V^E^^^l^^+^.-^^^+^J '(4-14) 

where R{uj) = {uj — J/2){uj — 3J/2){uj — 5J/2){uj — 3J). To recover what we should have 
obtained in the absence of renormalizations, we need only to replace R{uj) by {u — 2J){uj — 
4J)^(a; — 6 J) and {uj — 3 J/2) in the argument of Q hj {uj — 2J). 

We shall now solve the above equation numerically. Fig.|l^ shows our numerical results 
for the density of states for J = 0.5t on a cluster of 20 sites. Since our calculation is reliable 
only near the bottom of the spectrum and for relatively small J/t, one should note that 
the data on the high energy side of the spectrum does not make sense and is therefore 
omitted. In the following the discussion will be focused only on the lowest two peaks 
in the figure. First of all our numerical solution shows that the hole spectrum becomes 
dispersive with the introduction of vertex corrections and the ground state is located at 
k = (0,0). The density of states around the lowest peak therefore broadens in contrast 
to the (5-peak in the Shraiman-Siggia limit. Secondly, the peaks shift to lower energies 
when including the renormalization whereas the gap between the lowest two peaks becomes 



narrower as is shown in fig.|T2|a and can readily be explained from the form of R{uj). Also 
remarkable is the fact that the gap fits nicely a law ~ aJ°'^^ for < 0.4, with or without 
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the renormalizations (see the figure caption for the prefactor a). This is close to the J^/^- 
bahavior in the Shraiman-Siggia hmit. Let us now take a further look at the quasiparticle 
peak at the bottom of the spectrum. As was explained before, the nature of the quasiparticle 
band does not depend too much upon the renormalizations we have included as long as J/ 1 



is small compared to unity. In fig.|TpD we show quasiparticle residues for momenta at both 
the top and bottom of the band. With or without the renormalizations, the two residues 
are remarkably close to each other. The residue at the bottom of the quasiparticle band 
is always the smallest for > 0.3. Moreover, Zk=o shows a J^'^^-bahavior for Jz/t < 0.6. 
These results agree well with numerical findings by exact diagonalization0. In fact, the 
structure of the quasiparticle band can actually be obtained analytically now that only 
the leading vertex correction are included. With the dominant pole approximation^, we 
readily find that the quasiparticle band obeys: 

£k = ^0 — ^"^{(cosfca; + COS kyY + cos^k^ + ky) + cos(A;a; — ky)} , (4-15) 

where Eq is a constant. The quasiparticle is therefore seen to be confined within a narrow 
band of bandwidth J / and the minimum at A; = (0, 0). Finally, the ground state energy 
for the hole exhibits the J^'^^-behavior, simply following the string scenario. This should 
not be a surprise in view of the energy scale of the vertex correction near the bottom of the 
band. 



C. The coherent spectrum: spin fluctuations 

Finally we switch back on the spin fiuctuations which necessarily favor the delocalization 
of the hole. For simplicity, we shall take into account only the fiuctuations in the transverse 
spin channel. The interaction part in Hi is then composed of eq.( [4.4|) and the the Heisenberg 
exchange in the transverse direction. As a immediate consequence we now have to consider 
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also the process described by fig.p3 since the hole may easily face a "wrong-spin" neighbor 
created through spin fluctuations. It is certainly likely to hop to that neighbor in order 
to restore the local Neel ordering by absorbing the overturned spin. On the other hand, if 
the hole does hop to one of its Neel ordered neighbors, spin fluctuations will help to lower 
the energy by dispersing the newly created spinfiip through the medium. In summary, spin 
fluctuations facilitate the propagation of the hole either by creating a trap in front of it 
or by immediately eliminating the trace it leaves behind. In this section we shall limit 
ourselves only to this new mechanism, while neglecting Trugman's mechanism assuming 
that the latter plays only a relatively minor role. This neglect is in fact justified by the 
numerically small prefactor in front of the fc-dependent term in ( |4.15| ) {J/Az"^ = J/64). 



In order to understand the diagrammatic description, let us look at the case in which a 
downspin hole initially sits on the origin and is surrounded by a local Neel environment. As 
it hops to one of its neighbors, it turns to an upspin hole at the new position and creates a 
spinfiip at the origin. The flipped spin will now propagate in the medium because of spin 
fluctuations. Suppose that it approaches the hole again and ends up at one of its neighbors 
(on the even sublattice). For the reasons discussed above the hole would like to eliminate 
that overturned spin. This can happen in two different ways. It may either directly hop 
there, or it may hop one step more to another neighbor and let the newly created spinfiip 
be eliminated together with the previous one through spin fluctuations. In the first case, 
all it involves in the intermediate states is the propagation of a spinfiip from the origin to a 
site on the same even sublattice. The Feynman diagram for this process is therefore fig.p!3|a, 
where the propagation of the intermediate spinfiip is given by uj) from eq.( |3.1CI| ). The 

selfenergy of the hole is then determined by 

mk,uj + tO+)=t^Yl / -^Q^ik-q,uj-u;^ + t0+hl_^xuiq,^i) ■ (4-16) 

It differs from eq. ( [4.7D only by the fact that the intermediate particle-hole pair is now able 

24 



to propagate. On the other hand, in the second case we see an even number of particle-hole 
bubbles in the intermediate state between two sucessive hops of the hole. The process is 
therefore described by fig.p!3|b and we obtain: 



j:\{k,uj + tO+)=i^^J2j^g^k-q,u-Ui+t0^hk-qlkX2iiq,uj,) . (4.17) 

Some attention should be paid to the interaction vertex (namely 7fc) at the second hoping 
process. Its momentum dependence is determined by that of the hole in the final state {i.e. 
k). As for the dynamical susceptibility, it is given by the off-diagonal term in eg. ( p.lOD 
since we actually see the propagation of a pair of spinflips created by two successive hopings. 
An analogous analysis applies if the hole is initially surrounded by distorted spins. The 
corresponding Feynman diagrams are shown in fig.|TB|c and figjTB|d, from which we find: 

2z'^f^ r flu: 

mk, CO + 20+) = Yl J -q,u;-u;, + t0^hlx22{q, ^i) , (4.18) 



and 



Sj(fc, cu + zO+) = Yl / -q,u-uj^ + zO+)7fc7fc-<,Xi2(g, ^i) , (4.19) 



respectively. Finally we integrate over frequency and sum eqs. (|4.16| )- (H:.19|) to get the total 
selfenergy: 

Ei(fc,^) = |5^/(fc,q)6;|'(fc-g,a;-^,) , (4.20) 

Q 

f{k, q) = l-zHm^=L= + 1)7L, + {^=L= - Ihl - . (4.21) 

2 Vl-7q Vl-7q V^-Tq 

The three terms appearing in the dressed vertex function f{k, q) reflect the different 
scattering processes we have justed discussed. (The factor 2 in the last term comes from the 
degeneracy in the processes leading to and S*^). In particular f{k,q) shows unambigu- 
ously that for a hole on the "nesting line" 1^x1 + l^yl = tt (where 7fc = 0), the only nonzero 
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contribution comes from S". The hole having these particular momenta responds only to 
its Neel ordered neighbors. In addition, this hole only couples weakly to the dispersing low 
energy spin excitations, leading to f{k,q) ~ (q^ + qy)'^/\q\, which in turn translates into 
JmS ~ Lu'^ in two dimensionJii. In other words, quasiparticle behavior is well defined for 
these particular momenta. Secondly, quantum fluctuations play their full role in the hole 
dynamics only away from nesting wavevectors. Nonetheless, at these momenta none of the 
aforementioned scatterings by themselves support quasiparticle picture, as is clearly indi- 
cated by the individually divergent terms in (|4.21|) . Only from the sum of all the terms do 
we obtain f{k, q) ~ \q\ or JmS ~ uj'^ in the long wavelength limit and hence a quasiparticle 
behavior for the single hole. Finally, we notice that eqs. ([4.20| , |4.21|) are exactly the same 
as those obtained by Schmitt-Rink et al.i3 and by Kane et ahlli using the slave-fermion 
technique. Numerical calculations based on eqs.( [4.20D and (|4.21| ) have been carried out by 
several groupsii'0'§i and they have indeed observed a well-behaved quasiparticle peak near 
the bottom of the spectrum. In particular, the ground state is found to be at = (7r/2, 7r/2) 
where, as we now see, a quasiparticle picture is well estabilished. In ref. pO the numerical 
results obtained by the above groups were also compared with exact diagonalization results. 
Good agreement has been found concerning the k dependence of the quasiparticle residue 
as well as its variation with J. 

Before ending this section, we remind the reader that in studying the dynamics of a 
single hole we have assumed from the start the scattering processes involving energies of 
0{U) are irrelevant as far as low energy properties of the single hole are concerned. To 
see exactly how this comes about we can in principle follow the same reasoning as for half 
fllling. In the presence of spin fluctuations the selfconsistent RPA is then formed within the 
entire Hilbert space including scatterings involving double occupancies. The corresponding 



Feynman diagrams are those of flg. 13 except that the bare bubbles are replaced by those of 
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fig.|l]b. The final result is therefore again eg. ( |4.2CI|) and eq. ([4.21|) . In concluding this section, 
we also cite the recent calculation performed by Brenig and Kampflfl. These authors have 
studied the dynamics of a single hole in the presence of a SDW antiferromagnet which is 
then extrapolated to the large-U limit. When a selfconsistent RPA scheme is adopted (as 
is done here), they have observed, along with the incoherent background, a quasiparticle 
band with a significantly reduced bandwidth. Moreover, the quasiparticle residue is found 
maximal at k = (7r/2, 7r/2) and minimal at A; = (0, 0). These observations agree with those 
derived from eg. ( [4. 201) and eq. ([4.21|) . 



V. SUMMARY 

We summarize our main results in this section. We have shown that a transformation to 
a new basis facilitates the investigation of the crossover between weak and strong coupling 
limits of the Hubbard model. The new basis has been introduced via a Bogolyubov rotation. 
Its relation with the weak coupling method is seen from the fact that the transformation 
had previously been applied to the Hartree-Fock Hamiltonian which in turn served as the 
basis for further studies of fluctuation effects. Indeed, for the case of half fllling, we have 



explicitly demonstrated in section |T| the smooth crossover between two limits of Hubbard 
correlations. High energy scatterings have been shown to be irrelevant and are scaled away 
as U approaches the strong coupling regime. The correct spin wave mode we have derived 
in the flnal result is actuallly rather evident since the exact Heisenberg interaction in the 
transverse spin channel is contained in the low energy part of the effective Hamiltonian [i.e. 
eq.»). 



Within the new basis, we have further studied in detail the dynamics of a single hole 
in the antiferromagnetic background. Aside from the fact that we can diagrammatically 
distinguish the scatterings responsible for the low energy physics from those at high energy, 
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we are also able to separate in Feynman diagrams the incoherent from coherent physical 
processes. In the Ising limit, we have analyzed a series of vertex corrections responsible for 
the hole delocalization and have concluded that the leading vertex correction dominates as 
far as low energy physics is concerned. In particular, for small J/T quantitative agreement 
is reached with exact calculations on small clusters regarding the low energy physics, namely 
the ground state and the nature of the coherent motion for a single hole near the bottom 
of the spectrum. The transverse spin fluctuations further support a quasiparticle picture 
in the case of a single hole. We have derived a selfconsistent equation which recovers 
previous results for the t — J model which were obtained using slave-fermion methods. The 
numerical solution based on that selfconsistent equation also shows good agreement with 
exact solutions on small clusters. Finally, agreement is also found with results derived from 
direct extrapolation of weak coupling calculations recently done by Brenig and Kampfl 

An important ingredient in the present calculation is the two sublattice antiferromag- 
netic spin structure. This is certainly correct in the case of one hole. However, finite 
hole-doping will certainly destroy the antiferromagnetic spin structure on which our dis- 
cussions are based. The question concerning the stability of the quasiparticle picture and 
in particular its relevance to physical properties at finite hole densities therefore remains 
open. Further, the change of physical properties as one goes from weak to strong correla- 
tion at finite doping has also not been fully clarified. We hope to address this question in 
a forthcoming investigation. 
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APPENDIX A: 



In this appendix we explain the procedure which leads to eg. ( ^.2| ), i.e. the Hubbard 
Hamiltonian written within the reduced Brillouin zone (RBZ). 

By definition, the Fourier transform of a fermion operator at site j is expressed as 



^7 



rij 



where the summation over k covers the RBZ (i.e. {k^l + \ky\ < vr) and Q = (±7r,±7r). 
Ck+Q complements and together they run over every state in the whole Brilliuon zone 
once and only once. Correspondingly, the number operator is written as 

qk 

where again the summations for both q and k are confined to within the RBZ. In the mean- 
time, Ck+q and Ck+q+Q are complementary and cover the whole Brillouin zone centered at 
q. The above expression can then be written compactly as 

^^'e-*'^-j{^^(fc + g)^(fc) + + g)a^^(fc)} (A3) 

qk 

by defining a spinor ^E'J.(fc) = (C*fco-' ^fc+Qo-)- '^^^^ Hamiltonian can then readily be trans- 
formed noting that the following equality holds for the Fourier transform on sublattices: 

^Y!e-^'^- = 5,,o , (A4) 

r 

where the summation runs on either sublattice and q is defined within the RBZ. The 
?7-term in the Hubbard model can then be expanded as 

i kklq 

+ ^{(fc + q)(T^^|(fc)^|(fc'-q)(Tl^|(fc')} • (A5) 

Eq. (|2.2|) is thus derived. 
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APPENDIX B: 



In this appendix we detail the procedure which transforms the susceptibihty matrix, 
eq.(O), so that it is represented by the r-operators, i.e. in terms of Wannier functions. 
Following appendix A, the spin operator = Cj|Cjj at site j can be written within the 
RBZ as 



1 ^r^' 



(Bl) 



qk 



One can now Fourier transform the susceptibility matrix given in eg. within the RBZ 
and then write it in terms of r-operators by inverting eg. ( ^^.31 ). The final result for half 
filling is 
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As shown in eg . (|B3| ) , we need in general to calculate response functions with momentum 
dependent vertex functions fi{k) or f2{k) attached to either of the two external points. 



) 



(B4) 
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Nevertheless, for U ^ t, we can easily verify that x is dominated by the first term in 
eg. (|B3|) . The remaining terms are smaller by at least a factor 0{J/U). 
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FIGURES 

FIG. 1. (a) The RPA scheme for the spin-spin correlation function (xii) in the transverse 
spin channel. The dotted lines between bubbles denote the (72 interactions. The "bare" bubble 
here is the result of (b). (b) The interband multiple scatterings between a c-particle and v-hole 
which result in the renormalized "bare" bubble used in (a). The interaction vertex in the ladder 
are gi. 

FIG. 2. The 1st order nonzero contribution to X21 in eq.( |3.6| ). The "bare" bubble in the 
figure is the g'l -renormalized result given in fig.|l|b. As in Fig.|l] the dotted lines denotes the 52 
interactions. 

FIG. 3. The second order correction to the hole selfenergy. For a downspin hole, (a) is deter- 
mined by the first term in Hj, while (b) by the second term, (b) does not contribute in the absence 
of spin fluctuations (see text). The dashed line joining the bubble and the hole line denotes the 
0{t) interactions in Hj. 

FIG. 4. The selfconsistent 2nd-order correction to the selfenergy which accounts for the inco- 
herent background in the hole spectrum. 

FIG. 5. A typical selfenergy diagram in the Shraiman-Siggia limit. 

FIG. 6. The leading vertex correction in the Ising limit, (b) is an equivalent of (a), in which 
the particle-hole bubble is denoted by a dotted line; (c) The real space configuration corresponding 
to (a) or (b). 

FIG. 7. Nonzero vertex corrections at 10th order. The dotted line denotes a particle-hole 
bubble of opposite spins as in fig.^. Note that the spin of the c-particle in each bubble is the same 
as that of the intermediate third hole. 
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FIG. 8. Schematical representation in real space of the graphs in fig.^. 

FIG. 9. Maximally-crossed diagrams representing the coherent motion of a single hole in the 
Ising limit. FigjBjb and fig]^ are redrawn here in a different fashion to guide the eye. 

FIG. 10. The evolution of the energy in the intermediate states immediately after the nth 
hopping step: (a) for the leading vertex correction; (b) for the 10th order maximally crossing 
diagram (see fig.0a). 

FIG. 11. The density of states for a single hole in the Ising limit at J/t = 0.5, calculated on 
the 20~site cluster. A Lorentzian broadening with broadening parameter 5 = 0.01 is used, (a) 
With renormalzations as explained in the text, (b) without renormalizations. 

FIG. 12. (a) The gap A between the lowest two states as a finction of Jz (Ising limit), with 
renormalization (lower curve) and without renormalzation (upper curve) due to interactions in 
the longitudinal channel. The lines are fits to A = 2.610J^'^^ and A = S.ITSJ^'^^, respectively. 

(b) The quasiparticle residue for k = (0,0), denoted by ( + ) and k = (0,7r), denoted by ( * ), 
with renormalization (lower curves) and without renormalzation (upper curves). The lines are fits 
to Zfc=o = 0.511J°-S*^ and Zk=o = 0.662j0-86, respectively. 

FIG. 13. The selfenergy diagram in the presence of spin fluctuations: (a) the hole creates and 
then eliminates a spinflip at two successive hops; (b) it creates a spinflip at each of the two hops; 

(c) it first eliminates and then create a spinflip; (d) it eliminates a spinflip at each of the two hops. 
In (a) and (b) the hole faces a Neel surrounding in the initial state, while in (c) and (d) it sees 
one of its neighbors disturbed in the initial state. 
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